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OBSERVED  SOUND-SPEED  VARIABILITY  IN  THE  GULF  OF  CADIZ 

by 

James  Cairns  and  Lucio  Toma 


ABSTRACT 


f <• 


I 


Sound-speed,  fluctuations  have  been  examined  by  repeatedly  profiling 
temperature , conductivity , and  depth  over  long  tow  paths.  The  sound- 
speed  variability  due  to  internal  waves  conveniently  separates  from  that 
due  to  inhomogeneities , revealing  that  even  in  an  area  of  complex 
structure , internal  waves  are  the  legging  cause  of  sound-speed  fluctua- 
tions. At  very  low  wavenumbers  (<l(f  - aycle/m) , variability  due  to 
inhomogeneities  nearly  equalled  that  due  to  waves  only,  whereas  at 
higher  wavenumbers  they  contributed  only  about  one  fifth  as  much  vari- 
ability as  waves.  The  Munk  & Zachariasen  model  esttmate  of  wave-induced 
fluctuations  agrees  very  well  with  these  data,  and  allows  the  further 
computation  of  local  coherence  lengths.  A directional  dependence  of 
variability  is  found  near  the  coast.  Fluctuations  were  less  by  a factor 
of  five  looking  along  the  coast  than  looking  normal  to  it.  It  may  be  a 
general  result  that,  in  shoaling  areas,  sound-speed  fluctuations  measured 
parallel  to  the  shore  are  less  than  those  measured  at  right  angles  to 
the  shore. 
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x,z  Horizontal,  vertical  coordinates 

(")  Potential  quantities  referenced  to  sea's  surface 

Horizontal  wavenumber 
j Internal  wave  mode  number 

B Vertical  scale  length  = 1000  m 

oj.j  Inertial  frequency  - 0 05  cycle/h 

N(z),N  Buoyancy  frequency,  surface  extrapoiated  value 
< > Mean  value 
rms  ( ) Root  of  mean  squared  value 

S,0,p,P  Salinity  (%c).  Temperature  (°C),  Density  (g/cm3),  Pressure 
o = 1000  (p  - 1 0) 

a,b  Coefficients  of  thermal  expansion  and  haline  contraction 

a,6,y  = • Is  ’ C 
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INTRODUCTION 


Variability  of  the  oceanic  structure  on  vertical  scales  of  metres  to 
tens  of  metres  and  horizontal  scales  of  tens  of  metres  to  kilometres  is 
dominated  by  two  features  The  structure  is  characterized  by  inter- 
leaving globs  of  diverse  water  types,  which  sort  themselves  out  slowly 
under  the  influences  of  gravity,  friction,  and  diffusion  into  a system 
of  stable  layers.  They  are  evdent  as  the  steplike  irregularities 
typically  found  in  vertical  profiles  of  the  ocean's  temperature  and 
salinity.  Additional  variability  arises  from  the  continuous  wobbling 
and  straining  of  the  whole  layered  system  due  to  the  action  of  internal 
waves.  Both  waves  and  globs  contribute  to  irregularities  in  the  sound 
field,  and  although  their  relative  importance  is  not  clear,  internal 
waves  are  thought  to  be  the  dominant  source  of  sound-speed  fluctuations 
in  the  open  sea  [i ] . 

Differences  in  the  way  that  waves  and  globs  behave  allow  the  separation 
of  their  effects  They  are  thought  to  have,  first  of  all,  very  different 
time  scales,  globs  persisting  for  times  that  are  long  compared  with 
internal  wave  periods.  Secondly,  globs  are  characterized  by  a compensa- 
tion of  temperature  and  salinity  that  balances  out  horizontal  irregularities 
in  density  (but  not  in  sound  speed),  whereas  internal  waves  have  no  such 
compensation.  The  purpose  of  the  present  study  is  to  examine  experimentally 
the  separate  contributions  of  globs  and  waves  to  fluctuations  in  the 
sound  field 

Recently  an  experiment  was  staged  in  the  oceanographically  complex  Gulf 
of  Cadiz,  where  strong  currents  and  waters  of  widely  different  origins 
exist  in  close  proximity  (See  [2},  for  example).  One  would  therefore 
expect  this  to  be  an  area  where  sound  structure  anomalies  due  to  globs 
would  be  at  their  maximum.  In  spite  of  this,  the  sound-speed  variance 
in  the  Gulf  was  found  to  be  primarily  due  to  internal  waves. 


1 THE  EXPERIMENT 

The  experiment  was  conducted  from  SACLANTCEN's  R/V  MARIA  PAOLINA  G.  in 
June  1977,  Measurements  were  made  with  a towed  sensor  package  (TOB) 
developed  at  SACLANTCEN  [3].  The  TOB  instrumentation  consisted  of 
sensors  and  electronics  from  a Plessey  model  9400  CTD  modified  to 
include  a propeller-type  current  meter  in  a fourth  channel*.  It  was 
mechanically  designed  for  towing  at  modest  speeds.  A sketch  of  the  TOB 
and  its  towing  arrangement  is  shown  in  Fig.  1. 


Data  fanom  the.  cuAA&nt  meAzu  <vw  not  di&cuA6e.d  tn  tkei  pape-t. 
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While  towed  at  2 kn,  the  instrument  package  oscillated  in  depth  along  a 
taut  cable  between  120  m and  170  m.  Average  vertical  sensor  speed  was 
1 m/s;  instantaneous  speed  varied  somewhat  due  to  the  ship's  pitching. 
The  sensor  path  in  space  was  nominally  a 45°  zig-zag  line.  Only  the 
upward-travelling  profiles  were  analyzed,  giving  on  average  a profile  of 
50  m vertical  extent  every  100  m along  the  ship's  track.  The  routes 
along  which  the  tows  were  made  were  determined  by  satellite  navigation; 
they  are  shown  in  Fig.  2. 

Temperature,  conductivity,  and  pressure  were  simultaneously  sampled  8 
times  per  second  and  recorded  as  16  bit  binary  coded  numbers  on  magnetic 
tape.  Their  least  count  precisions  are  ±0.002SC,  ±0.002  ms/cm*,  and 
±0.034  m (equivalent  water  depth),  respectively.  Time  constants  for  the 
temperature,  conductivity,  and  pressure  probes  are  0.4,  0.1,  and  0.1  s 
respecti vely . The  data  were  filtered  to  minimize  the  effects  of  the 
pitching  ship,  and  time  shifted  to  account  for  the  diverse  sensor  time 
constants  (see  Appendix  A). 


1 .1  Data 

The  data  set  is  derived  from  about  350  km  of  tow  and  has  been  examined 
primarily  for  its  sound-speed  characteristics;  a more  detailed  discussion 
of  the  physical  structure  of  the  water  and  its  dynamics  will  be  given  in 
a separate  paper  [4] . 

It  is  assumed  throughout  the  discussion  that  the  structure  can  be  considered 
frozen  in  time,  so  that  at  a 2 kn  tow  speed  only  spatial  changes  are 
seen.  This  is  thought  to  be  a good  assumption  for  the  internal  wave 
field  [5]\  it  should  be  an  even  bettrr  assumption  for  the  glob  field. 

A profile  in  the  upper  kilometre  (Fig.  3)  gives  an  idea  of  the  gross 
sound-speed  structure.  It  was  obtained  from  a deep  CTD  cast  made  at  the 
end  of  tow  segment  5. 


The  oceanic  environment  encountered  along  the  tows  ranged  from  very 
quiet  areas,  where  both  internal  wave  and  glob  activities  were  low,  to 
regions  of  extreme  structural  convolutions  with  internal-wave  levels 
typical  of  the  open  sea.  Figure  4 shows  a section  of  profiles  demon- 
strating a quiet  region  along  track  8.  Each  profile  on  the  figure  is  a 
graph  of  one  of  the  variables  versus  depth;  successive  profiles  have 
been  plotted  by  shifting  the  scale  to  the  right  by  a fixed  distance.  On 
each  profile  the  depths  of  encounter  of  particular  values  of  the  variables 
were  registered.  Lines  passing  from  profile  to  profile  through  these 
depths  give  the  changing  elevation  of  the  isopleths  along  the  tow. 


Figure  4 shows  few  inversions  of  the  variables  w'th  depth,  and  the  up 
and  down  oscillations  of  the  iso-potential  density  surfaces  (i.e.  the 
internal  waves)  are  small. 


★ 

Tkz  ^.ntz'incuU.o  naZ  unZt,  AtzmznA  (4),  zqiM.vaZznt  to  thz  psizvZoaiZy- 
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FIG.  2 AREA  OF  EXPERIMENT,  SHOWING  SHIP'S  TRACKS 
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FIG.  3 SOUND- SPEED  PROFILE  DERIVED  FROM  A CTD  CAST  TAKEN  AT  THE 
END  OF  TOW  SEGMENT  5 
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In  contrast,  Fig.  5,  from  tow  segment  4,  depicts  one  of  the  most  active 
regions.  The  fields  of  temperature,  salinity,  and  sound  speed  appear  to 
have  been  folded  into  complex  patterns.  Inversions  in  these  fields  were 
tens  of  metres  high  and  several  kilometres  in  horizontal  extent.  The 
potential  density  field  on  the  other  hand,  which  was  presumably  affected 
only  by  wave  activity,  showed  no  evidence  of  these  distortions,  although 
the  wave  activity  itself  was  higher  here  than  in  the  section  shown  in 
Fig.  4.  Obviously  in  this  area  there  was  considerable  variation  of  the 
sound  speed  due  not  directly  to  internal  waves  but  to  globs. 


2 SEPARATION  OF  EFFECTS 

As  an  aid  to  understanding  sound-speed  variability,  the  contribution  due 
to  internal  waves  was  separated  from  that  due  to  globs,  using  a scheme 
worked  out  in  collaboration  with  Walter  Munk  of  the  Scripps  Institution 
of  Oceanography.  The  following  discussion  is  limited  to  fluctuations 
along  a horizontal  path,  as  this  is  the  only  situation  for  which  there 
were  observations.  We  make  the  assumption  that  all  the  observed  depth 
fluctuations  of  surfaces  of  constant  potential  density  are  due  to  internal 
waves. 

It  is  useful  to  consider  three  cases. 


Case  1 — waves  only,  no  globs 

In  the  presence  of  an  internal  wave  field  with  displacement  z, , a particle 
momentarily  at  reference  depth  zQ  has  come  from  a rest  position  zQ  - ?. 

The  resulting  sound  speed  at  zQ  is  estimated  to  be 

C(x,zQ)  = C0  + £ c(x),  [Eq.  1] 

where  CQ  is  the  average  sound  speed  along  z = zQ.  Effects  due  to  wave- 

induced  particle  velocities  and  pressure  changes  have  been  ignored,  an 
approximation  shown  to  be  feasible  by  Munk  & Zachariasen  [l] . 

To  estimate  the  wave-induced  variability  from  the  data,  where  both  waves 
and  globs  may  be  present,  we  write 

C(x,zQ)  = c0  + M 4(x) , 


with  (dC/dz)  the  local,  average,  vertical  sound-speed  gradient.  For  the 
wave-induced  sound-speed  fluctuation  along  zQ,  this  gives 


<5Cw(x’z0)  = ({£)  (x>  e(x). 


[Eq.  2] 
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Case  2 — globs  only,  no  waves 


In  the  absence  of  motion,  variability  along  a lateral  path  is  due  entirely 
to  globs.  Their  contribution  to  the  fluctuating  sound  speed  in  this 
case  would  be 


<5Cg(x,Z0)  - C(x,z0)  - C0.  [Eq.  3] 

An  equivalent  measure  of  6Cg  can  be  made  even  in  the  presence  of  waves, 

however,  if  observations  are  made  along  a wave's  surface.  Measurements 
along  the  wavy  surface,  when  adiabatic  effects  are  removed,  are  equiva- 
lent to  measurements  made  along  its  rest  depth  in  the  absence  of  mction 
A surface  of  constant  potential  density  moves  with  the  waves.  Making 
observations  along  the  potential  density  surface  with  rest  depth  z 
gives  glob-induced  sound-speed  fluctuations  of 


<5Cg(x ,zq)  = C(x,  c(x) ) - CQ, 


[Eq.  4] 


/N 

where  C = C - yz  . 
o o o 


The  relative  importance  of  salinity  and  temperature  to  glob-induced 
sound-speed  fluctuations  is  seen  as  follows.  By  definition,  along  an 
iso-potential  density  surface 


a60  = b6S , 

where  in  our  case  a»  0.2  x 10~3(°C)_1  , and  bw0.8  x 10~3(%o)~*. 
Furthermore  we  may  write 

6Cg  = 6©( ac/ae)  + 6S(3C/3S), 


The  ratio  of  saline-induced  to  thermally-induced  sound-speed  fluctua- 
tions along  the  isopycnal  is  then 


. 6S(ac/as)  = as 
6e( ac/ae)  ba 


[Eq.  5] 


where,  for  this  case,  a = C"‘(3C/36)  = 2.5  x 10'3(CC)-1  and 
6 = C_1(aC/3S)  = 7.7  x 10~H(%o)~1,  giving  c - 0.08. 
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Case  3 - waves  and  globs 

This  is  the  total  sound-speed  variability  one  observes  along  z = z in 
the  presence  of  both  waves  and  globs,  which  is  simply 

6C(x,z0)  = C(x,zQ)  - CQ  [Eq.  6] 


3 ANALYSIS  AND  RESULTS 

A 

The  data  series  c(x),  C[x,  c(x)],  and  6C(x,zo),  allowed  examination  of 

the  records  for  total  sound-speed  variability  and  for  the  separate 
contributions  due  to  waves  and  globs.  The  c(x)  series  was  obtained  by 
finding  the  elevation  of  a particular  iso-potential  density  surface  on 
every  upward  movement  of  the  TOB,  or  roughly  every  100  m along  the 
ship's  track.  These  series  were  then  resampled  at  an  even  100  m sampling 
length  using  linear  interpolation.  At  each  resulting  point  (x,  s(x)), 

C was  determined*,  yielding  the  second  series  of  interest.  Finally 
C(x,zQ)  was  found  at  the  same  horizontal  distances  along  the  constant 

depth  z = zQ.  Using  these  variables  determined  from  the  measurements, 

and  Eqs.  2,  4,  and  6,  the  additional  series  6Cw(x,zo),  6Cg(x,zQ),  and 

6C(x,zq)  were  formed  for  each  segment  of  the  tow  pattern  shown  in  Fig.  2. 

Examples  of  the  data  series  are  given  in  Fig.  6.  They  are  from  tow 
segment  3,  which  appears  qualitatively  to  be  in  an  area  with  a typical 

level  of  activity.  These  series  were  then  multiplied  by  C to  give 
series  of  fractional  sound-speed  fluctuations. 


3.1  The  Spectra 

The  series  are  not  stationary  for  spectral  computations,  so  they  were 
pre-whitened  and  the  spectra  later  adjusted  in  order  to  compute  the 
curves  shown  in  Fig.  7**. 

We  see  that  at  all  wavelengths,  the  total  fractional  variability  (6C/C) 
was  approximately  equal  to  (6C.7C).  At  wavelengths  less  than  5 km, 

variability  due  to  waves  exceeded  that  due  to  globs  by  roughly  a factor 
of  2.  Variability  due  to  globs  decreased  somewhat  faster  with  wave 
number  than  did  that  due  to  waves,  but  they  both  varied  approximately  as 

_ 2 

k„  . This  was  expected  for  (6C ,/C),  as  demonstrated  in  the  next  section. 

I ! W 


★ ^ 

Both  C and  C were  computed  by  uiing  equation*  &rom  [6],  When  calcula- 

A 

ting  C,  z it  * at  equal  to  zero  and  temperature  t&  replaced  by  potential 
temperature . 

irk 

The  high  wavenumber  portion * o&  the  observed  spectra  are  *omeiwhat 
Autpect  due  to  vertical  filtering  Uee  App.  A) . 
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MODEL  ESTIMATED  SPECTRUM  IS  FROM  MUNK  & ZACHARIASEN  11] 
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3.2  Model  Comparison 

Munk  and  Zachariasen  [i]  find  for  the  power  spectrum  of  fractional 
sound-speed  fluctuations  due  to  internal  waves* 


F(6C  /C)(kH 
w 


where 


.» ■<(£)) 


6CwA  2n~  h/No)B’  kH  jH(j) 


[(kH)2  + (WB"  <yN0)2]5 


. 1 T-* 


h(j)  = (j2  + j*r  //_,  (j2  + j*)'  > 

j=i 

£ h(j)  • i. 
j=i 

£ (j2  + j*  k i. 

j=i 

Then  summing  over  all  inodes  yields  the  wavenumber  spectrum 


[Eq.  7] 


[Eq.  8] 


/ (5  C 2v  ^ 

F(6C  /C)(kH^  \ ("r)/27T  Z 

w ' / 4-1 


kj  jH(j ) 


j=l  C(kH)2+  (%jB"  ui/N0)2]2 


[Eq.  9] 


where  B = 1000  m,  j*  = 3,  and  we  use  our  local  value  of  wi 
(=  0.05  cycle/h) . 

From  Munk  and  Zachariasen  we  make  the  following  estimates: 
N = NQe'2/B;  (N  = 3 cycle/h),  and 
rms(6Cw/C)  = rms(6Cw/C)0(N/N0)3/2 


[Eq.  10] 


OuA  notation  <u>  itightly  dt^eaent  &Aom  theAd.  We  aie  cycttcal 
aatheA  than  c -OicuJtan.  tmvenumbeAi  &oa  eaie  o£  compaction  u)tth  the 
obieAved  data. 
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wi  th 


rms(6Cw/C)0  = (3.0  x 10'6)  uN*  rmsU^g'1 ; 
rms(c0)  = 7.3  ni  from  [5]  and 


U = (a/a)  s(Tu). 


[Eq.  11] 


[Eq.  12] 


In  Eq.  12,  s = (1  + cTu),  where  c is  as  defined  previously  (Eq.  5),  and 
where  the  Turner  number  is 


Tu  = 3 1 (a  3z^  [Eq.  13] 

CiQ  _ o _ , 

Our  gradients  are  typically  = -2.2  x 10  {%0)m  , and 

||=  -1.4  x 10"2(°C)m"  \ yielding  Tu  = 0.63,  and  s = 2.8. 

From  Eqs . 12,  11,  10,  we  calculate  the  model  estimates  y = 22.0; 

rms(6C  / C ) = 4.5  x 10  4 ; and  rms(6C  /C)  = 3.6  x 10’4  The  estimated 

WO  w 

mean  squared  value  from  the  model  is  then 


’•3* ,0'’ 


Replacing  this  in  Eq.  9 gives  us  the  Munk  & Zachariasen  model  spectrum 
of  internal-wave-induced  fractional  sound-speed  fluctuations  shown  in 
Fig.  7.  It  is  to  be  compared  to  the  spectrum  of  observed  wave-induced 
fluctuations,  with  which  it  agrees  remarkably  well. 


3.3  Coherence  Scales 

To  the  extent  that  the  Munk  & Zachariasen  model  represents  the  Gulf  of 
Cadiz  environment,  we  may  further  use  it  to  estimate  the  local  hori- 
zontal (Lh)  and  vertical  (L y)  coherence  scales.  The  model  gives 


ff(N0/o)i  )B 

8j*[log  (N/coj  )-l/2] 


= 2.2  km. 


BNq/N 


= 0.14  km. 


[Eq.  14] 


[Eq.  15] 


r 


These  scales  are  z-dependent  throughout  N(z).  The  effect  of  y obs 
will  make  the'  actual  coherence  distances  somewhat  less  than  these, 
and  the  above  values  should  probably  be  viewed  as  upper  limits. 


3.4  The  Nearshore  Region 

Spectra  of  sound-speed  fluctuations  along  the  individual  tow  segments 

were  all  more  or  less  alike,  with  one  exception.  Measurements  parallel 

to  the  bottom  contours  in  the  shallow  (500  m)  region  traversed  by 

track  8 (Fig.  2)  were  markedly  quieter  than  the  others.  An  example  of 

the  track  8 data  is  shown  in  Fig  4.  The  spectrum  of  fluctuations 

from  this  path  is  compared  in  Fig.  8 with  the  overall  spectrum 

from  all  paths.  We  see  that  sound-speed  fluctuations  along  segment  8 

are  about  a factor  of  5 lower  across  nearly  the  whole  spectrum,  reflecting 

the  fact  that  both  glob  and  internal  wave  activities  were  low  along  the 

path.  This  is  thought  to  be  due  primarily  to  anisotropy  of  the  wave 

field  in  this  shoaling  water  region,  with  wave  crests  aligned  more 

parallel  to  the  coast  than  perpendicular  to  it.  If  so,  it  is  likely 

to  be  a general  result  that  in  the  nearshore  region  sound-speed  fluctuations 

parallel  to  the  coast  will  be  less  than  fluctuations  along  paths  normal 

to  the  coast. 


CONCLUSIONS 


Sound-speed  fluctuations  have  been  examined  by  repeatedly  profiling 
temperature,  conductivity,  and  depth  over  long  tow  paths.  The  sound- 
speed  variability  due  to  internal  waves  conveniently  separates  from  that 
due  to  inhomogeneities,  revealing  that  even  in  an  area  of  complex  structure, 
internal  waves  are  the  leading  cause  of  sound-speed  fluctuations.  At 

very  low  wavenumbers  (<10~ “cycle/m  ),  variability  due  to  globs  nearly 
equaled  that  due  to  waves  only,  whereas  at  higher  wavenumbers  they 
contributed  only  about  one  fifth  as  much  variability  as  waves. 

The  Munk  & Zachariasen  model  estimate  of  wave-induced  fluctuations 
agrees  very  well  with  these  data,  and  allows  the  further  computation  of 
local  coherence  lengths. 

We  find  a directional  dependence  of  variability  near  the  coast.  Fluctua- 
tions were  less  by  a factor  of  five  looking  along  the  coast  than  looking 
normal  to  it.  We  think  it  may  be  a general  result  that  in  shoaling 
areas  sound-speed  fluctuations  measured  parallel  to  the  shore  are  less 
than  those  measured  at  right  angles  to  the  shore. 
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FIG.  8 SPECTRUM  OF  FRACTIONAL  SOUND  SPEED  FLUCTUATIONS  ALONG  TOW  TRACK  8 
(shallow  water  alongshore  path)  COMPARED  WITH  OVERALL  SPECTRUM 
FROM  ALL  PATHS.  DEGREES  OF  FREEDOM  ON  OVERALL  SPECTRUM  RANGE 
FROM  12  to  96,  AND  ON  SEGMENT  8 SPECTRUM  FROM  8 to  30 
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APPENDIX  A 


There  are  some  inherent  problems  in  making  measurements  from  a pitching 
ship  using  sensors  wUh  finite  and  diverse  time  responses.  These  are 
especially  bothersome  when  information  from  multiple  sensors  must  be 
combined  to  calculate  another  variable  such  as  density  or  sound  speed. 
Pingree  [A.l]  gives  a good  description  of  the  situation  A sensor 
lowered  from  an  equilibrated  position  requires  some  time  to  notice  it 
has  moved  If  lowered  uniformly  through  a smooth  gradient  it  eventually 
parallels  the  true  trace  but  lags  behind  it.  If  its  descent  slows  or 
stops,  however,  it  will  continue  to  equilibrate  and  its  trace  will  curve 
over  to  join  the  true  one-  Alternate  slowing  and  speeding  of  the  descent 
results  in  a wobbly  trace  that  lags  the  true  smooth  trace;  the  effect  is 
reversed  on  the  ascent.  In  principle  one  can  remove  these  effects  if 
the  sensor  has  sufficiently  well  known  and  well  behaved  characteristics. 
But  in  practice  their  removal  is  quite  difficult,  especially  as  they 
effect,  say,  computed  density. 

Wiggles  were  found  in  our  density  profiles  with  a regularity  of  about 
once  each  six  seconds  (or  each  six  metres  vertically).  The  power 
spectrum  of  vertical  sensor  speed  derived  from  the  first  difference 
record  of  the  pressure  sensor  was  also  sharply  peaked  at  1/6  Hz,  indica- 
ting the  pitching  response  of  the  ship  to  the  encounter  of  the  ocean 
waves 

We  have  minimized  these  effects  in  our  data  at  the  expense  of  vertical 
resolution  by  filtering  all  of  the  data  with  a filter  whose  vertical 
high  wavenumber  cutoff  is  0.07  cycle/min  (14  m vertical  scale).  To 
compensate  for  the  phase  shift  induced  by  its  longer  constant,  the 
temperature  record  is  then  shifted  back  in  time  by  0.34  s («1  time 
constant)  relative  to  conductivity  and  depth  This  time  shift  was 
determined  empirically  by  minimizing  the  small-scale  variance  in  com- 
puted density  profiles.  All  the  variables  including  the  calculated 
density  and  sound  speed  are  then  very  well  behaved,  but  we  are  con- 
strained to  look  at  features  with  large  vertical  scales. 


